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1. Introduction 

In the problem of bounding the size of codes in compact homogeneous 
spaces, Delsarte's polynomial method gives rise to the most powerful uni- 
versal bounds on codes. Many overviews of the method exist in the lit- 
erature; see for instance Levenshtein [1]. In this note, which extends our 
previous work [2] we develop a functional perspective of this method and 
give some examples. We also discuss another version of the functional ap- 
proach, a linear algebraic method for the construction of polynomials for 
Delsarte's problem. Our main results are new constructions of Levenshtein's 
polynomials. 

Let X be a compact metric space with distance function r whose isom- 
etry group G acts transitively on it. The zonal polynomials associated with 
this action give rise to a family of orthogonal polynomials V{X) = {P^} 
where k = 0, 1, . . . is the total degree. These polynomials are univariate 
if G acts on X doubly transitively (the well-known examples include the 
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Hamming and Johnson graphs, their (/-analogs and other Q-polynomial 
distance-regular graphs; the sphere 5"^^^ G M'') and are multivariate other- 
wise. 

First consider the univariate case. Then for any given value of the degree 
K = i the family V{X) contains only one polynomial of degree i, denoted 
below by Pi. Suppose that the distance on X is measured in such a way 
that t{x,x) = 1 and the diameter of 3C equals —1 (to accomplish this, a 
change of variable is made in the natural distance function on X) . We refer 
to the model case of X = S"*^^ although the arguments below apply to 
all spaces X with the above properties. Let {f,g} ^ J^i fgd/J, be the inner 
product in L2{[—1, where dfj.{x) is a distribution on [—1, 1] induced 

by a G-invariant measure on X. Let ^ (•, 1) be the moment functional 
with respect to dfi. We assume that this distribution is normalized, i.e., 
that ^(1) = 1. 

Let C be a code, i.e., a finite collection of points in X. By Dclsartc's 
theorem, the size of the code C whose distances take values in [— l,s] is 
bounded above by 

|C|< inf /(l)//o, (1) 

where 

^ = ^s)^{f:f{x)<0,xe[-l,s]; /o > 0, > 0, i = 1, 2, . . . } (2) 

is the cone of positive semidefinite functions that are nonpositive on [—1, s] 
(here /j = {f, Pi) /{Pi, Pi) are the Fourier coefficients of /). 

2. Functional approach 

The choice of polynomials for problem (l)-(2) was studied extensively in 
the works of Lcvenshtein [4-6] . In this section we give a new construction 
of his polynomials and their simplified versions. 

2.1. Notation. 

Let V be the space of real square-integrable functions on [—1, 1] and let Vk 

bo the space of polynomials of degree k or less. Let pi = Pi/ {Pi, Pi), i = 
0,1;... be the normalized polynomials. The polynomials {pi} satisfy a 
three-term recurrence of the form 



xpt = ajpj+i -I- biPi + ai_ipi_i, 

i = l,2,.. . = 0,po = 1; a_i = 0. 



(3) 
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In other words, the matrix of the operator x : V ^ V (multiplication by the 
argument) in the orthonormal basis is a semi-infinite symmetric tridiagonal 
matrix, called the Jacobi matrix. Let Xk = Eko x where E/. = pioiy^y^ , 
and let Jfe be the (fc + 1) x (fc + 1) submatrix of J, 



Jk 



60 ao ... ■ 

oo 61 ai ... 
ai 62 02 ... 

dk-i 

. flfe-i h . 



Example 2.1. (a) For instance, let X be the binary n-dimensional Ham- 
ming space. Then Pi{x) = fcj("/2(l — x)), where ki{z) is the normalized 
Krawtchouk polynomial. The polynomials Pi{x) are orthogonal on the finite 
set of points {xj = 1 — {2j/n),j = 0, 1, . . . ,n} with weight w{xj) = (p2~" 
and have unit norm. In this case. 



a, = {l/n)^{n - + 1), b, = 0, < i < n. (4) 

(6) Let X be the unit sphere in d dimensions. Then Pi{x) are the normalized 
Gegenbauer polynomials; in this case 



y (n + 2t){n + 2i — 2) 



It is well known [7, p. 243] that for A; > 1 the spectrum of Xk coincides 

with the set X^+i = {xk+1,1, ■ ■ ■ ,Xk+i,k+i} of zeros of Pk+i- Below we 
denote the largest of these zeros by Xk+i- Let 

k 

Kk{x,s) = ^Pi{s)pi{x) (5) 

be the fc-th reproducing kernel. By the Christoffel-Darboux formula, 

(a; - s)Kkix,s) = ak{pk+i{x)Pk{s) - Pk+i{s)Pk{x)). (6) 

In particular, if s G "Xk+i then XkKk{x, s) = sKk{x, s). Note that Kk{x, y) 
acts on Vfc as a delta-function at y: 

{Kki;y)J{-)) = f{y). (7) 
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2.2. Construction of polynomials. 

Without loss of generality let us assume that /(I) = 1. Then (1) is equiva- 
lent to the problem 

sup{nf),/e*}- 

Let us restrict the class of functions to Vn- By the Markov-Lucacs theorem 
[8, Thm. 6.4], a polynomial f{x) that is nonpositive on [—1, s] can be written 
in the form 

fn{x) = {x-s)g^-{x + l)4il or fn{x) = {x + l){x - s)g'^ - (f)l 

according as its degree n = 2fc + 1 or 2A; + 2 is odd or even. Here g, (pi G 
Vk,(f>2 G Vfc+1 are some polynomials. Below the negative terms will be 
discarded. We use a generic notation c for multiplicative constants chosen 
to fulfill the condition /(I) = 1. 

2.2.1. The MRRW polynomial. 

Restricting our attention to odd degrees n = 2fc + 1, let us seek f{x) in the 
form (x— s)g^. Let us write the Taylor expansion of in the "neighborhood" 
of g. Let h €Vk he a. function that satisfies \\h\\ < e for a small positive e 
and the condition h{l) = 0. We obtain 

J^{{x - s){g + hf) = T{{x - s).g2) + {{x - s){g + h),g + h)- {{x - s)g,g) 

= T{f) + T'(h) + y2{r'h,h), 

where J^' = 2(x — s)g,T" — 2{x — ,s) are the Frechet derivatives of T. 
This relation shows that for / to be a stationary point of J^, the function 
g should satisfy dJ^ = 2{g, [x — s)h) = for any function h with the above 
properties. First assume that s = Xk+i- Then by (6), a stationary point of 
is given hy g = Kh{x, s), and we obtain / in the form 

fn{x) = C{X - s){Kk{x,s))'^. 

Since /o = 0, conditions (2) are not satisfied; however, it is easy to check 
that they are satisfied ii Xk < s < Xk+i- For all such s, the polynomial /„ 
is a valid choice for problem (1), yielding 

\C\ < , , Kl{l, s). (8) 

' ' - akPk+i{s)Pk{s) ' ^ ' 

The polynomial /„ was used by McEliece el al. [9] and Kabatiansky and 
Levenshtein [10] to derive their well known upper bounds on codes. 
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2.2.2. Levenshtein polynomials, n = 2fc + 1. 

So far in our optimization we did not use the condition h{l) = 0. To use it, 
let us write h = {1 — x)hi,h\ G Vfc-i and repeat the above calculation. We 
find that stationary points of should satisfy 

dJ^^-^ = 2{{x - s)g, (1 - x)hi) = 0, 

where J^^~^ {■) = J ■ {1 — x)diJ, is the moment functional with respect to the 
distribution dfj,^~\x) = (1 — x)d^{x). A stationary point of J^^~^ is given 
by the reproducing kernel K^{x, s) with respect to this distribution: 

k 

=^pr(s)pr(a;), (9) 

where {p~{x), i ~ 0, 1, ... } is the corresponding orthonormal system. To 
find the polynomials (x) observe that 

:F^-\p-p-) = Hv-{x)p-{x){i - x)) = 

is satisfied for p^ (x) = Ki{l, x) / {aiPi+i{l)pi{l)Y^^ . Indeed, if j < z then 
the function (1 — x)Ki{l, x) is in the subspace spanned by pi+i,pi and thus 
is orthogonal to Kj{l,x). To conclude, the function sought can be taken in 
the form 

f-{x)=c{x-s){Ki;{x,s)f. 

2.2.3. Levenshtein polynomials, n = 2k + 2. 

In this case wc seek the polynomial in the form /„ = [x — s){x + 1),9^. The 
necessary condition for the stationary point takes the form J^^[{x — s)gh) = 
{{x — s)(l — x'^)g,h) = 0. Prom this, g = K^{x,s) where the kernel 
is taken with respect to the distribution diJL^^\x) = (1 + a;)(l — x)dii{x). 
The corresponding orthogonal polynomials pf{x) are also easily found: up 
to normalization they arc equal 

pf{x) = Ki{x, - Ki{x, 

Then 

.ftix)=cix-s)ix + l)iK^ix,s)r. 

Let x'^ (x^) be the largest root oip^{x) (resp. ofp^(x)). Then G $ 

if x^ < s < x^^i and fij,+2i^) e * if x^+i <s< x^^^. 
Remarks. 

1 . The polynomials /~ , were constructed and applied to coding the- 
ory by Levenshtein [4-6]. Polynomials closely related to them were studied 
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in a more general context in the works of M. G. Krein et at: see Krein and 
Nudelman [8]. The orthogonal systems {p^},{pf} are sometimes called 
adjacent polynomials of the original system {pi}- 

2. The stationary points found above are not true extremums because 
the second differential of the functionals J- . J-'^^\ J-^^^ is indefinite: for 
instance, (fij-{g) = 2{(x — s)h,h). Nevertheless, the polynomials fn^fn 
have been proved [11] to be optimal in the following sense: for any n > 1 
and all / e deg f <n 

Hfn) > Hf)- 

3. Asymptotic bounds derived from (1) relying upon the polynomials 
fn,fn^fn Coincide. For the finite values of the parameters, better bounds 
are obtained from ,fn- 

3. Spectral method 

This section is devoted to a different way of constructing polynomials for 
the Delsarte problem. The ideas discussed below originate in the work of 
C. Bachoc [12]. They were elaborated upon in an earlier work of the authors 
[2]. 

We develop the remark made after (6), namely that for any i > 1, 
Kk{x, Xk+i,i) is an eigenhmction of the Jacobi operator X]^. Since if/c(x, s) 
is a good choice for the polynomial in Delsarte's problem, it is possible to 
construct polynomials as eigenvectors of Xk as opposed to the analytic ar- 
guments discussed above. In particular, Kh{x,s) arises as an eigenfunction 
of the operator = Tfc(s) defined by 

Tk -.Vk ^ Vk 

(f) ^ Xk(j) + Pk4>kPk 

where pk = akPk+i{s) / Pk{s) ■ Indeed, using (5) and (6) we obtain 

{Tk - s)Kk{x, s) = {Xk - s)Kk{x, s) + akPk+i{s)pk{x) = 0. 

On account of earlier arguments we should choose the polynomial for prob- 
lem (1) in the form F{x) = {x — s)p{x) where f{x) = f{x, s) is an eigen- 
function of Tk- The positive definiteness condition of / can be proved using 
the Perron-Frobenius theorem; for this we must take / to be the eigenfunc- 
tion that corresponds to the largest eigenvalue of Tk . This condition defines 
the range of code distances s in which the method is applicable. 

A variant of this calculation was performed in [2] to which we refer for 
details. The difference between [2] and the argument above is that there we 
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took pk = afePfe+i(l)/pfc(l). This has the advantage of defining indepen- 
dently of s but leads to a bound of the form 



which is generally somewhat weaker than (8). Using the function F defined 

above we can improve this to recover the estimate (8). 

We note that this argument does not depend on the choice of the func- 
tional space; in particular, the kernels K^,K^ arise if the operator 
is written with respect to the basis of the corresponding adjacent poly- 
nomials ({p^} or {pf}) and their generating distribution. To conclude, 
Levenshtein's polynomials and bounds on codes can be derived within the 
framework of the spectral method. 

Example 2.2. Consider again Example 2.1(a). The adjacent polynomials 
up to a constant factor that does not depend on i are given by [5, p. 81] 



where fc, {z) for instance denotes the degree-i normalized Krawtchouk 
polynomial associated with the (n — l)-dimensional Hamming space. The 
Jacobi matrix for the basis pj can be computed from (4) as follows. 
Since 



we find that the coefficients of three-term recurrence for the family {p^ } 
are 

a, = (l/n)V(n-fc - 1)(A: + 1), 6^ = -1/n, i = 0, 1, . . . . 

Constructing the operator as described above, we obtain K^{x, s) as its 
eigenfunction. A similar construction can be pursued for the function K"^. 

The approach outlined above has two advantages. First, it enables one to 
obtain simple estimates of the largest eigenvalue of X/. which is important 
in verifying the condition f{x) < 0,x G [— l,s]. The second advantage 
is a more substantial one: this method can be extended to the case of 
multivariate zonal polynomials when the analytic alternative is not readily 
available. This case arises when the space 3i is homogeneous but not 2-point 
homogeneous. Worked examples include the real Grassmann manifold Gk,n 
( [12]; the Pi are given by the generalized /c-variate Jacobi polynomials) and 



\c\< 



4afcPfc+i(l)pfe(l) 



(10) 



p-{x) = H"-'\z), pf{x) = U"-'\z) for z = ^(1 - x) - 1 
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the so-called ordered Hamming space [3] . We provide a few more details on 
the latter case in order to illustrate the general method. 

Let Q be a finite alphabet of size q. Consider the set Q'"'" of vectors 
of dimension rn over Q. A vector x will be written as a concatenation of 
n blocks of length r each, x = {xn, . . . , Xir] . . . ; a;„i, . . . , Xnr}- For a given 
vector a; let Cj, i = 1, . . . , r be the number of r- blocks of x whose rightmost 
nonzero entry is in the ith position counting from the beginning of the 
block. The r- vector e = (ei, . . . , e^) will be called the shape of x. A shape 
vector e = (ei, . . . , e^) defines a partition of a number N < n into a sum 
of r parts. Let eo = n — e^. Let Ar,„ = {e e (Z+ U {0})'" : J2i — 
be the set of all such partitions. The zonal polynomials associated to Q""'" 
are r-variate polynomials -P/(e), /, e G A^.n of degree n = /j. They are 
orthogonal on Ar,„ according to the following inner product 

^ Pf{e)P,{e)w{e)=0 {f ^ g). 

The weight in this relation is given by the multinomial probability distri- 
bution 

w{ei,...,er) = n\'[l-^ {Pi = q'~''~'^{q - l),i = 1, . . . ,r;po = 5"''), 
i=o 

so the polynomials Pf{e) form a particular case of r-variate Krawtchouk 
polynomials. 

Let X G Q*"'" be a vector of shape e. Define a norm on Q*"'" by setting 
w(a;) = zcj and let dr{x, y) = w{x — y) be the ordered Hamming metric 
(known also as the Niederreiter-Rosenbloom-Tsfasman metric). We note 
that in the multivariate case there is no direct link between the variables 
and the metric. For instance, for the space Q*"'" the polynomials (as well 
as relations in the corresponding association scheme) are naturally indexed 
by shape vectors e while the weight is some function e. 

The Delsarte theorem in this case takes the following form: The size of 
an {n, M, d) code C C Q'"-" is hounded above by M < inf^£$ ./(0)//o, where 

* = {fix) = f{xi, . . . ,Xr) =fo + Yl /-^«(^) ■■fo>0,fe>0{ej^ 0); 

/(e) < Ve s.t. ^ ie, < d] 

i=l 

The argument for the univariate case given in this section can be repeated 
once we establish a three-term relation for the polynomials Pf{e). Let be 
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the column vector of the normalized polynomials Pj ordered lexicographi- 
cally with respect to all / that satisfy J2i fi = and let F{e) be a suitably 
chosen linear polynomial. Then 

F(e)P«(e) = A«P.+i(e) + B^V^{e) + Al_^¥^^i{e) 

where A^^.B^, are matrices of order ('^^I^^) x (""""^^p^) and s = 1,0, re- 
spectively. The elements of these matrices can be computed explicitly from 
combinatorial considerations. This gives an explicit form of the operator 
5k = i?K o F{e) in the orthonormal basis. Relying on this, it is possible to 
derive a bound on codes in the NRT space of the form (10) and perform ex- 
plicit calculations, both in the case of finite parameters and for asymptotics. 
The full details of the calculations are given in [3] . 
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